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We consider holographic cosmological models of dark energy in which the infrared cutoff is 
set by the Hubble's radius. We show that any interacting dark energy model with a matter 
like term able to alleviate the coincidence problem (i.e., with a positive interaction term, 
regardless of its detailed form) can be recast as a noninteracting model in which the holo- 
graphic parameter c 2 evolves slowly with time. Two specific cases are analyzed. First, the 
interacting model presented in [l[ is considered, and its corresponding noninteracting ver- 
sion found. Then, a new noninteracting model, with a specific c 2 (z) expression, is proposed 
and analyzed along with its corresponding interacting version. We constrain the parame- 
ters of both models using observational data, and show that they can be told apart at the 
perturbative level. 



I. INTRODUCTION 

Nowadays the Universe appears to be undergoing a phase of accelerated expansion [i}. This can be explained 
through modified gravity theories [S| and by general relativity (GR) modulo introducing an exotic component, dubbed 
dark energy (DE). We only know for certain that DE is endowed with a hugely negative pressure (of the order of 
its energy density) and that it should be distributed rather evenly across space; see [j to learn about the state of 
the art. Considering a spatially flat Friedman-Robertson- Walker (FRW) metric in the frame of GR the Friedmann 
equations read 



(pM + px) , (1) 
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M~ 2 

H = ^-(pM + px + P), (2) 

where Mp = (87rG)~ 1,/2 is the reduced Planck mass, P = wpx, and px and pM indicate the DE and dark matter 
(DM) energy densities, respectively. 

Thus far, the most successful candidate for DE is the cosmological constant A, which together with cold dark 
matter and radiation form the standard cosmological model, ACDM. Regrettably this model suffers from two main 
problems: (i) the cosmological constant value, A, is many orders of magnitude lower than that expected for the 
vacuum energy density from quantum field theory, and (ii) the coincidence problem. 

The latter can be alleviated by introducing a suitable interaction between DE and DM, which appears consistent 
with observations from relaxed galaxy clusters @43- For interacting models, the conservation equations are 

pM + 3HpM = Q and p x + 3H(1 + w)p x = —Q ■ (3) 
where Q is the interaction term. 

Whatever the nature of DE it seems reasonable that it respects the holographic principle. The latter asserts that 
the number of relevant degrees of freedom of a system dominated by gravity must vary as the area of the surface 
bounding the system @|. In addition, the energy density of any given region should be bounded by that ascribed to a 
Schwarzschild black hole that fills the same volume [9(. Mathematically this condition reads px < Mp L~ 2 , where L 
stands for the size of the considered region (i.e., the infrared (IR) cutoff). This expression is most frequently written 
in its saturated form 

3Mpc 2 ... 
= • ( 4 ) 

Here c 2 is a dimensionless parameter -very often, assumed constant- that summarizes the uncertainties of the theory 
(such as the number particle species and so on); the factor 3 was introduced just for mathematical convenience. 

When dealing with holographic DE one must first specify the IR cutoff. With the lack of clear guidance different 
expressions have been adopted. The most relevant ones are the Hubble radius, -ff -1 , see e.g. [l(J , the future 
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events horizon Rh = J°° * that suffers from a severe circularity problem, and the Ricci length, Rcc, i.e., 

L = Rcc = {H + 2H 2 )~ 1/2 -see e.g. [12h16| ]. The rationale behind the latter is that it corresponds to the size of the 
maximal perturbation leading to the formation of a black hole [TtI ] . 

It has been argued that an IR cutoff defined by the Hubble radius, L — H~ , can not lead to an accelerated 
Universe, since DE evolves as H 2 , that is, as pressureless matter, therefore producing a decelerated expansion [ill ]. 
However, if DM and DE interact via Eqs. @, with Q > 0, the present epoch of accelerated expansion can be realized 

Differentiating the expression px = 3Mj>c 2 H 2 , that describes DE density for Hubble holographic models, and 
using Eqs. ([2]) and (02), the equation of state (EoS) parameter 

(5) 



3(1 - c 2 )p x H 

follows. Note that Q must be positive, otherwise w would be positive, and the Universe would not accelerate. 
Further, if Q were negative, the second principle of thermodynamics would be violated [lJj|. A comprehensive study 
of holographic DE models can be found in |20t | - 

An alternative approach to simultaneously solve the coincidence problem and describe the late time acceleration 
was taken in [2^]. To get an evolving energy density ratio, r = pm/px, with the Hubble scale as IR cutoff, the 
c 2 parameter was promoted to vary with time, though slowly, in such a way that (c 2 )' > 0. In that scenario, the 
holographic bound is progressively saturated in the sense that c 2 ^.^ = constant > 0. 

Likewise, it was shown that a variable c 2 significantly alleviates the coincidence problem also for cosmological 
model with nonvanishing spatial curvature [22I ]. An extended analysis of nonsaturated holographic DE models with 
nonvanishing spatial curvature was recently presented in [23|. There, it was shown that a variable c 2 is compatible 
with current precision data, which favor a small Q^t Du t high enough to have significant cosmological consequences. 

Moreover, a nonsaturated holographic bound provides a transition from a decelerated to an accelerated era not 
only for standard holographic DE models, but also for a class of generalized holographic DE models in which the 
gravitational coupling is promoted to a time dependent quantity [241 ]. 

In [25| the role played by a variable c 2 parameter in a model with the IR cutoff set by the Hubble's length 
was discussed in connection with the behavior of the effective EoS of time variable cosmological constant models. 
A parametrized expression for c 2 was proposed and tested using observational data. Notice, however, that this 
parametrization of c 2 , leads to a negative DM density in the future. 

In 126] a modified holographic DE model with a variable c 2 was considered in which the IR cutoff was set by the 
Ricci scale, L = Rcc- Three different parametrizations were proposed and tested with cosmological data leading to 
consistent results. 

Recently, a comprehensive analysis of variable c 2 holographic DE models for the three widely used infrared cutoff 
scales, namely the Hubble's length, the Ricci's length and the particle horizon length, was presented in [2?| . It was 
argued that the c 2 term appearing in the conventional formula for holographic DE should not be assumed constant 
in general, except for the particular case of DE with the Ricci's length as infrared cutoff. 

In what follows, quantities referring to holographic models with variable c 2 will be noted by a tilde. They obey 

p M = 3Mp(1-c 2 )H 2 and px = 3Mj>c 2 H 2 . (6) 

as well as 

Pm = -3Hp~M and px = —3H(1 + w)p~x , (7) 
i.e., we assume that their energy densities conserve separately. 

In this paper, we study some general features of holographic DE models where the IR cutoff is defined by the Hubble 
radius. By considering both points of view, we shall show that, in general, identical cosmological backgrounds can be 
described by an interacting holographic DE model where the holographic parameter c 2 is a constant or, alternatively, 
by a noninteracting holographic DE model in which the c 2 parameter depends weakly on time. In spite of the global 
evolution in both scenarios being the same, the energy densities, the EoS parameters and so on, can behave rather 
differently. 

We remark that in any holographic interacting model defined at the Hubble scale, independently of the nature 
of the interaction, DE and DM share the same dependence on H, and thus present the same background evolution. 
By rewriting Eqs. ([3]) as 

pu + 3H ( 1 - — ^— \ p M = and p x + 3 H (l + w+-§—) p x =0 (8) 

and using Eq.([5} both, DM and DE, are seen to share the same effective EoS parameter, w e // = ~ zjnr~ ■ We wm 
return to this later on. 

We shall also r-consider the coincidence problem from the viewpoint of noninteracting models where the energy 
densities ratio, f = pM /px, is not a constant. 
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The plan of the paper is as follows. In Sec. [TTJ we show the equivalence at the background level between interacting 
holographic DE models and noninteracting holographic DE models in which the holographic parameter c 2 depends 
on time. In Sec. Mil we briefly introduce the cosmological tools (e.g. luminosity distance, angular distance, and 
so on) used to characterize the models studied in this work. In Sec. IIV1 a specific model is considered, namely, 
the interacting DE model of Ref. Jl,]. Its corresponding non-interacting version is derived and their Lagrangian 
formulation is presented. In Sec. [V] another model, in which a novel expression for c 2 is proposed, aimed to 
alternatively describe the late behavior of both DE and DM energy densities, is presented and contrasted with the 
concordance ACDM model. In section ["VTl we show that the background equivalence is broken at the perturbative 
level. Finally Sec. I VIII summarizes our main results. 



II. BACKGROUND EQUIVALENCE BETWEEN INTERACTING AND c 2 MODELS 

To show that every interacting model can be considered as a noninteracting one (with a c 2 parameter varying in 
time) at the background level (i.e., both interpretations share the same Hubble function), we must first verify that 
both DE and DM energy densities are positive. By Eqs. ([6]), this condition implies < c 2 < 1. 

From Eqs. H), @ and (02) we obtain 

H = ~^H 2 (l + c 2 w) . (9) 

For the noninteracting model, from Eq. (01) we have that pm = pMod remains always positive. 

We must also verify that px is positive, i.e., that c 2 > 0. Since by hypothesis interacting and c 2 models share 
the same H(z), Eq. ((9J implies 

^P = ^4, (io) 

where z is the redshift (z — - — 1). Thus, if the noninteracting DE has a negative EoS parameter w(z), since 
w(z) < and c 2 > 0, then px and c 2 (z) will always be positive. 

Since pm > and c 2 < 1, by differentiating Eq. (J6j2) , using Eqs. ([5]), (0) and (fT0|) . and recalling that Q > 0, we 
obtain 

(£ 2 )'=i^Q>0. (ll) 

PM 

As shown in [27|. this condition must be fulfilled for thermodynamical reasons. Thus, we conclude that the Hub- 
ble function of any holographic interacting model with Q > and c 2 constant also corresponds to a holographic 
noninteracting model (Q — 0) of DE (w < 0), with c 2 obeying < c 2 (t) < 1 and (c 2 )' > 0, and vice versa. It is 
also noteworthy that despite both models being equivalent at the background level, they share the same H(z) and 
their energy components evolve diversely. In the interacting case, the Hubble function will never have c 2 as a free 
parameter, but it will be multiplied by the constants in the interacting term Q, as can be seen by using Eq. (0 
in Eq. (0). This means that, as long as 0.x = c 2 , the fix o parameter cannot be fitted because we have neither 
observational nor theoretical constraints on the interaction, and so the dependence on them is degenerated. On the 
other hand, in the c 2 model it has the fixed value, tlx o = <?(z — 0). Consequently, at background level it is not 
possible, in principle, discriminate interacting models from c 2 ones. Nevertheless, as we shall see in Sec. IVI1 they 
are distinguishable at the perturbative level. 



III. OBSERVATIONAL CONSTRAINTS 



To constrain the free parameters of the Hubble holographic models presented below (Sec. IIVI and [V)) . we shall 
use observational data from SN la Union2 set (557 data points) [28(], BA O 1291 ] . the acoustic scale I a [30|], gas mass 
fractions in galaxy clusters as inferred from x-ray data (42 data points) [31| and the Hubble rate (15 data points) 
[32H35I]. As usual, the likelihood function is defined by C oc exp(— x 2 /2). The best fit parameter values can be found 
by minimizing the sum xLai = Xsjv + Xbao + Xi A + xl-ray, + XHubble- 

A. SN la 

We contrast the theoretical distance modulus 

Pth(zi) = 51og 10 (^^j + Mo , (12) 
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whereito = 42.38 — 51og 10 h, with the observed distance modulus ^i bs{zi) of the 557 SN la compiled in the Union2 
set [28J|. Here h = ^ with Ho in km/s/Mpc. The latter assembly is much richer than previous SN la compilations 
and has some other advantages, especially the refitting of all light curves with the SALT2 fitter and an enhanced 
control of systematic errors. In ()12[) di = (1 + z) J* E ^'- p ) denotes the Hubble-free luminosity distance, with p the 
model parameters. To eliminate the effect of the nuisance parameter ^to we resort to the method of [36| . 



B. BAO 

Pressure waves originated from cosmological perturbations in the primeval baryon-photon plasma produced acous- 
tic oscillations in the baryonic fluid. These oscillations have been unveiled by a clear peak in the large scale correlation 
function measured from the luminous red galaxies sample of the Sloan Digital Sky Survey (SDSS) at z = 0.35 [37| 
as well as in the Two Degree Field Galaxy Redshift Survey (2dFGRS) at z = 0.2 [29J]. These peaks can be traced to 
expanding spherical waves of baryonic perturbations with a characteristic distance scale 

D v (z B Ao) = 

wheredA(z) = (1+z)' 1 f* ^C. Data from SDSS and 2dFGRS observations yield D v (0.35) /D v (0.2) = 1.736±0.065, 
a nearly model independent value [2!|. The most constraining values r s {z. k )/D v {zi), for Zi — 0.2, 0.35, 0.278, obtained 
in pffll l3Sl]. are not model independent but obtained just for fiat and open ACDM and 10CDM models. 
The sound horizon radius is given by 

and the sound speed before decoupling is c 3 — ^1 + ^t 2 -) ' ', where Qb and fi 7 are the fractional density 
parameters of baryons and photons, respectively. 



zbao 



H{z B ao 



■ <1a{zbaoY 



(13) 



C. Acoustic scale Ia 

The cosmic microwave background (CMB) power spectrum is sensitive to the distance to the decoupling epoch 
via the peaks location. More specifically, with CMB data one can measure two distance ratios. The first one is the 
acoustic scale, I a, described by the angular distance to the decoupling surface divided by the sound horizon radius 
at that time 

; dA{z*) . . 

lA = 7T ^)> (15) 

where <1a(z) and r s (z) are comoving quantities. 

The other one is the CMB shift parameter, R, defined by the ratio of the distance to decoupling to the Hubble 
horizon, H~ 1 (z ir ), at that time 

R = %/£Wffo / -jt^ . (16) 
Jo H(z) 

The latter expression assumes a negligible DE density at decoupling, and should be modified when used in models 
other than ACDM or uiCDM [30L [33. |40|] . For this reason (our interacting model has a constant DE-DM ratio), we 
will just use the acoustic scale , I a, to constrain the models. However, to describe the Universe before decoupling 
one should include also the baryon and radiation energy component (see Appendix yQ . 



D. Gas mass fraction 



As is well known, a very useful indicator of the overall cosmic ratio Qbaryona/^M , nearly independent of redshift, 
is the fraction of baryons in galaxy clusters, f ga a\ see [4l]]. This quantity can be determined from the x-ray flux 
originated in hot clouds of baryons, and it is related to the cosmological parameters by f gaa tx d 3 / 2 , with d,A the 
angular diameter distance to the cluster. 

We used measurements by the Chandra satellite of 42 dynamically relaxed galaxy clusters in the redshift interval 
0.05 < z < 1.1 [3l|. In fitting the data we resorted to the empirical formula 

KAjbjz) Q B0 / 4 CDM V /2 n7) 
l9as[z) - 1 + s(z) Qmo \ d A ) [ ' 
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(see Eq. (3) in Ref. [3l[]) in which the ACDM model serves as a fiducial model. Here, the parameters K, A, 7, b(z) 
and s(z) model the abundance of gas in the clusters. We set these parameters to their respective best fit values for 
the considered model, but in the intervals defined in [3lj . 



E. History of the Hubble function 

Recently, high precision measurements by Riess et al. at z — 0, from the observation of 240 Cepheid variables of 
rather similar periods and metallicities [33 ] . as well as measurements by Gaztanaga et al, at z = 0.24, 0.34, and 0.43 
[33I ]. who used the BAO peak position as a standard ruler in the radial direction, have somewhat improved our 
knowledge about H(z). However, at redshifts above, say, 0.5 this function remains largely undetermined. Yet, to 
constrain the holographic model we have considered these four data alongside 11 noisier data in the redshift interval 
0-1 z < 1.8, from Simon et al. [13] and Stern et al. [3^], obtained from the differential ages of passive-evolving 
galaxies and archival data. 



IV. MODEL 1 

Let us now consider the holographic interacting model studied in [l|] to construct its equivalent c 2 (t) model. In 
the former the IR cutoff is also set by the Hubble length and the interaction term was taken as Q = SAHopM, 
with A a semipositive definite constant, related to the constant decay rate, F, of DE into DM by A = ■ The 
corresponding Hubble function is 

H = H (a + (1-A)(1 + *)*) , (18) 

and the equation of state parameter 

w = -^li' (19) 

follows from Eq. ((5| and the fact that Q.x = c 2 and Q.m = 1 — c 2 . Notice that here the density parameters, 
fli = 3h j'l H 2 (i = M, X) are constant. The model fits reasonably well the observational data and is consistent with 

the age of the old quasar APM 08279+5255 [i3 |. something that ACDM is not at 1<t confidence level. 
This model presents an unexpected similarity to the Chaplygin gas model [43^ . To see this, bear in mind that in any 
Hubble holographic interacting model DM and DE share the same effective EoS parameter, w e ff = —Q/(3Hpm). 
Multiplying it by the total energy density we get the total pressure 

P = w eff p = -3?AMpHopi , (20) 

with A a semi-positive definite constant. This expression is formally identical to the pressure of the generalized 
Chaplygin gas, P = ~/3p~ a [ij|. As is well known, for a < and ft > 0, it may imply instabilities since the squared 
adiabatic sound speed {c? sa = -?) is negative . Becaue of the interaction, we can take account of non adiabatic 
processes, and so consider an effective speed of sound. In Appendix [B] we find a Lagrangian formulation with a 
standard scalar field 4> f° r both models (the interacting and the <?) with an effective speed of sound given by c 2 = 1. 
Notice that while in the interacting case r = pmjpx is a constant, in the noninteracting one, f may vary with 
time. We expand H 2 (z) and assume that the term proportional to (1 + z) 3 , corresponds to the usual matter term in 
Friedmann's equation, and identify the remainder as the DE energy density. Thus, 



M~ M~ 3 

-^p M = (l-^) 2 (l + z) 3 and rp x = A 2 + 2A(l-A)(l + z)?, (21) 

and from Eqs. (J6j2) , (JT5J and J2TJ it follows that 

_ 2 = 2A{l-A)(l + z)i+A 2 ^ 
A)(l + 2)i) 2 

In consistency with the findings of 27] , c 2 never decreases and tends from below to a constant value in the far future, 

c!Loo=0, c 2 ^ = 1 - (1 - A) 2 , c 2 ^_x = l, (23) 

where < A < 1 from observations. Table U shows the best fit values of the parameters and their la errors. Table 
[H] presents the x 2 of model 1 and ACDM, obtained by fitting each data set independently, and the total x 2 - The 
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obtained \ 2 P er degree of freedom (dof) is \ 2 /dof — 1.00. Notice that (1 — A) 2 = Qmo an d its value (ss 0.17) is about 
6<7 lower than the one reported by Komatsu et al. pjpl ]. However, this value is reached by using the last scattering 
sound horizon as a standard ruler; i.e., it is not observed directly but by integrating the background evolution. So 
truly, just a global background evolution is obtained, which in the case of ACDM, gives the value of Qmo mentioned 
above, but in other models, as in this one, it can vary. Figure [TJ shows the la and 2a confidence regions and the best 
fit value of the free parameters of the models. 



Model 


^10 


Ho 


A 


Interacting holographic 


0.73 


69.4 ±1.7 


0.588 ±0.004 


c 2 holographic 


0.830 ±0.003 


69.4 ±1.7 


0.588 ±0.004 


ACDM 


0.720 ±0.003 


7i-5±l:l 





Table I. Values of the parameters of the models obtained by constraining them with the observational data described in Sec. 
IHII fix o is not a free parameter in any of the two holographic scenarios, but it is included for the sake of comparison with the 
conventional ACDM model. Notice that despite the fact that holographic models are described by the same Hubble function, 
they have different values of £lxo- This point does not discard any model, since the value from WMAP7 |30| is obtained using 
the last scattering surface as a standard ruler, i.e., integrating the ACDM background evolution from the last scattering surface 
to the present time. Thus, if the Hubble function is the same for different models, geometric data alone can not distinguish 
among them. The Ho values are given in km/ s/Mpc. 



Model 


XsN 


Xbao 


xl 


X.x — ray 


X 2 H 


Xtot 


Xtot/dof 


Model 1 


554.8 


1.7 


0.3 


44.9 


11.3 


613.0 


1.00 


ACDM 


542.7 


1.2 


0.7 


42.3 


8.8 


595.7 


0.97 



Table II. x 2 values for the holographic model of Sec. IIVI and the ACDM model. Each of them has two free parameters [A 
and Hq the holographic, and Qxo an< i Hq the ACDM). 
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69 - 

s 

X 
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67 - 



66 1 ' 1 ' 1 

0.58 0.585 0.59 0.595 

A 



Figure 1. l<r and 2cr confidence regions for the parameters A and Ho of model 1. The plot holds for both scenarios, the 
interacting and the c 2 (t) (noninteracting). The dot indicates the best fit values. 



Although H(z) coincides with the corresponding expression in PJ, pM, px, and w{z) do not. The functional form of 
w coincides with that of the EoS parameter for the interacting case - cf. Equation (1191) -. 

^4 Mo 

w = -nx-lT (24) 

but has a different time dependence since Six 7^ fix, as the left panel of Fig. [2] shows. Notice that in the interacting 
case, w crosses the phantom divide line (w = —1). However, the w a ff defined in the line below Eq. ([8|, does not 
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cross the said line. The pressure, obtained by multiplying Eq. (|24|l by Eq. ((6]2) and using Eq. (|18[1 . described here 
for later purposes, reads 

P = w~p x = -ZAMlHl [A+ (1 - A){\ + z)i^ . (25) 

As shown in the right panel of Fig. [2] the coincidence problem is solved (i.e., r is constant) in the interacting case 
(solid green line). By contrast, in the c 2 model (thin dot-dashed red lines), it is not solved but is much less severe 
than in the ACDM model (thick short dashed blue line). 




Figure 2. Left panel: EoS parameter for the interacting (to thin line, and w e ff thick line), the c 2 and the ACDM models. 
Right panel: energy densities ratios, r = Pm/px, versus 1 + 2 for the ACDM, the interacting and the c 2 models. All the graphs 
were plotted using the best fit values of the parameters, shown in table [T] Solid (green) lines are used for the interacting case, 
thin dot-dashed (red) lines for the c 2 model, and thick short dashed (blue) for ACDM. The lu region of the parameters is also 
plotted, but due to the very small errors, it results are nearly inappreciable. 



V. MODEL 2 

In this model, DM and DE evolve separately (i.e., Q = 0) but the holographic parameter c 2 varies slowly with 
time. To have < c 2 < 1, and (c 2 ) ' > 0, we use the parametrization 

(26) 



l + f (l + z) e 



where fo = EMfi an(1 e is a semipositive definite constant (notice that c 2 (z = 0) = flxo)- In this case, the Hubble 
function 



H = Ho\/«Mo(l + z) 3 + 0*0(1 + z) 3 " £ (27) 

coincides with the Hubble function of the spatially flat wCDM model with w = — |. It obviously reduces to the 
ACDM model for e — 3. If we consider the Hubble function in Eq.(|27p as coming from an interaction between DE 
and DM, by Eq. §5$ the interacting term would be 

Q = -3c 2 wp M H , (28) 

where the EoS parameter of the interacting case is w = — . a ,, , , — 777 . 

The variation of c 2 breaks the holographic dependence of DE density. But as for the rate of variation of c 2 we must 
have (c 2 ) '/c 2 = efo(l + z) e _ff/(l + fo(l-|-2:) £ ) ^ H, it may be considered as a small variation in the level of saturation. 
For the last inequality to be fulfilled, f(z) = ef2.Mo(l + zY/(Qxo + Qaio(1 + z) e ) < 1. This was not always true in 
the past, but the maximum value of f(z) , that is monotonously decreasing, is f(z — > oo) = e. Using the best fit 
value for e, shown in Table [TTT1 its maximum variation rate is of the order of the expansion rate, and the model can 
still be considered holographic. The \ 2 values obtained by fitting the model with the different data sets, are shown 
in Table [TTTl The \ 2 P er dof obtained with the fit of all the data sets together is x 2 /dof — 0.97. 
The left panel of Fig. shows that the EoS parameter w (thin solid green line) crosses the phantom divide line 
but the effective one Weff (thick solid green line) of the interacting case does not. The noninteracting case, as we 
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have mentioned before, is just a 10CDM model with w = — § (thin dot-dashed red line). The right panel of Fig. 
[3] shows that the interacting model (solid green line) solves the coincidence problem, and in the case of the c 2 one 
(thin dot-dashed red line), it overlaps the ACDM line (thick short dashed blue line), since iv — 1. Figure U depicts 
the la and 2a regions for the parameters fo and Ho (left panel) and fo and wq (right panel). Both panels are the 
same for both scenarios, the interacting and the c 2 , however, fo is only related to the DE and DM densities in the 
c 2 description; in the interacting case it has no physical meaning. In both cases today's value of the EoS is the same 
wq — wo — — | , as the left panel of Fig. [3] shows. 




. i i 

1 + z 



Figure 3. Left panel: EoS parameter for the interacting (w with thin line, and w e ff with thick line), the c 2 and ACDM 
models. Right panel: energy densities ratios, r = Pm/px, versus 1 + z for the ACDM , the interacting and the c 2 models. 
Notice that the energy densities ratio of the c 2 and ACDM models practically overlap. All graphs were drawn using the best 
fit values of the respective parameters, shown in Table Hill Solid (green) lines are used for the interacting case, thin dot-dashed 
(red) lines for the c 2 , and thick short dashed blue for the ACDM. The la region of the parameters is also plotted. 



Model 




Ho 


fo 


e 


Interacting holographic 


0.73 


71.5 ±2.6 


0.370 ±0.013 


2 97+ - 16 

z - y ' -0.14 


c 2 holographic 


0.730 ±0.007 


71.5 ±2.6 


0.370 ±0.013 


9 qy+0.16 
z - y ' -0.14 


ACDM 


0.720 ±0.003 









Table III. Best fit values of the free parameters of the models. In the c 2 holographic scenario, Clxo is obtained from the free 
parameter fo, so the model has only three free parameters. In the interacting one, for fixo we use the value obtained in [30ll . 
The Hq values are given in km/s/Mpc. 



Model 


XsN 


Xbao 


xf A 


2 2 

X.x — ray Xh 


Xtot 


Xtot/dof 


Model 2 


542.7 


1.2 


0.0 


41.5 8.8 


594.2 


0.97 


ACDM 


542.7 


1.2 


0.7 


42.3 8.8 


595.7 


0.97 



Table IV. x 2 values for the holographic model studied in Sec. |V]and for ACDM model. In the former, the free parameters 
are fo, e and Hq. In the latter the free parameters are two, Qxo an£ I ^0- 



VI. EVOLUTION OF THE SUBHORIZON PERTURBATIONS 

We have seen that in model 1 (Sec. IIV|) . the interacting version is pre ferred over the noninteracting one because 
&mo ~ 0.17, too small as compared with results of Komatsu et al. [30(]. In model 2 (Sec. [V]), both seem to be 
compatible with observations at the background level. A further study should be done, in general, to establish which 
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0.37 0.38 



.39 0.4 



0.35 0.36 0.37 0.38 0.39 0. 



Figure 4. Left panel: The la and 2<r confidence regions for the parameters fo and Hq of model 2. Right panel: The same 
regions for the parameters wo (equivalent to wo and derived from e = — 3wo) and fo. Notice that fo is the ratio of the energy 
densities just for the c?(t) (noninteracting) case; for the interacting case it has no straightforward physical meaning. Central 
dots indicate the best fit values. 



scenario is observationally favored. Here, we make a simple perturbation study, just considering matter perturbations 
inside the horizon. This study just intends to illustrate that a perturbative analysis can permit us to distinguish the 
interacting from the noninteracting scenario, despite the fact that they share the same H{z). 

Once the Universe becomes matter dominated, the anisotropic stress due to neutrinos will be negligible and, in 
the Newtonian gauge, it will suffice just one single gravitational potential -say <f>- to determine the flat metric element 

ds 2 = -(1 + 2(j})dt 2 + a 2 (l - 2(j>)dx i dx i . (29) 

In the interacting case, the energy-momentum tensors of DM and DE are not independently conserved. For the 
matter component, by perturbing the conservation equation T M U = Q M , where Q M = Qu M (with the four 
velocity of the matter component, and no momentum transfer in the DM reference frame is assumed) 0, [4j|, the 
equations of motion for the energy density and the velocity divergence [9 — ik 3 Vj) perturbations, with the speed of 
sound, c 2 M — 0, and the EoS parameter, wm = 0, are 

S M = -—+H + — + — {fl>-SM), (30) 
a pM Pm 

k 2 

9 M = -H9 M + —<t> , (31) 

and for the DE component T£ u = Q^-, where = —Qu M . Since DE can be considered a scalar field, the effective 
speed of sound c 2 x = 1, the dynamical equations for perturbations are 

5 X = - (1 + w) — - m (1 - w) 5 X - 3aH [3H (l - w 2 ) + wx] ^§ (32) 



+ 3(l + w)<j>- -2- 
px 



<t>-5 x + SaH (1 - w) ^ 



9 X = 2H9 X + tz 1 r — Sx -—<f>- 7— (fljwr - 29 x ) , (33) 

(1 + w) a a (1 + w) p x 

while the Fourier transformed time-time and space-time components of the Einstein equations are 

'„2 



fc-' 



■3H(4> + H<t>) = -4-KG^2 Pi 6i, (34) 

i 

— + H4>) = i-KG^2( Pi + Pi) 9, . (35) 

i 

A more detailed derivation of the perturbation equations can be found in [45j. As we are concerned with subhorizon 
scales, we just consider the case in which k 3> aH, and so, the Newtonian limit of Poisson's equation (|34p is just 

k 2 

—(j) = ~4ivG (pmSm + pxSx) ■ (36) 
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From Eq. (|36[> and bearing in mind that every single energy component obeys ^f^-pi < H 2 , one sees that the 
gravitational potential (and its derivatives) can be neglected when compared with the density perturbations. After 
all this, Eqs. (O, (EH, (E3) and JMI simplify to 



VM 

a 



(37) 



M = —HOm + —<f> , (38) 

5x = - (1 + w) — - 3H (1 - w) 8x + — {QSx - SQ) , (39) 
a px 

Ox = tt-^-Sx - jr-^y- (&M - 29 x ) , (40) 
(l+w) a (1 + w) p x 

Notice that in the interacting version of model 1, SQ = 3AHqpmSm , but in model 2, since Q includes a dependence in 
w and H, things are more involved. However, we assume, as in [4^], that the product 3c 2 wH is just an approximation 
to a time (but not position) dependent interaction rate, so there are no perturbations on it, and SQ — —3c 2 wH pmSm ■ 
In the c 2 scenarios of both models Q — and SQ = 0. 



A. Initial conditions 



To solve numerically these four coupled differential equations, we must choose some initial conditions for the 
density and velocity perturbations. We set them at z% = 1000. We impose the potential to be a constant (we have 
seen that it and its time derivative are much smaller than density perturbations), and using the perturbed Einstein 
equations, (|35p and (|36|l . we find the density and velocity initial conditions 

fc = — , (pMi5M(Zi,k) + pxiSx(zi,k)} , (41) 

^ \ -L -r Zi) 

k 2 (p= ctn 3 .„ (pMi0 M (zi,k) + pxi0x{zi,k)) . (42) 
2(1 + Zi)H 

for each model. Since we obtain, numerically, that the evolution of perturbations is nearly independent of the 
wavenumber k in the range 0.001 hMpc~ x < k < 0.1 hMpc -1 , which includes all the interesting scales un- 
der consideration, we shall assume that initially DE perturbations are proportional to the DM perturbation, i.e., 
Sx{zi,k) = a Sivr(zi, k) and 9x(zi,k) — /30M(zi,k), with a a nonnegative constant. 



1. Interacting version of model 1 

We have freedom to normalize the matter density contrast as 8m (zi) = 7^7-, and then by Eq. (|35[) and (|36[) . 

we find 6m {zi) = — n^°"n^ ) • ^ or ^ ne com P onen t we take different options for both a and /3. The 

dashed (green) lines of the left panel of Fig. [5] correspond to the noninteracting scenarios. The solid (blue) line 
is for the ACDM model, shown for comparison purposes. Each different dashed (green) line depicts various initial 
conditions, from top to bottom, a = and /3 = 0, a = 1 and /3 = 1, and a — 1 and f3 = 10. Notice that the 
interacting scenarios exhibit a dependence on the initial conditions chosen for the DE component. The evolution of 
perturbations depend on the initial conditions on Ox but not so much on Sx- 



2. c 2 version of model 1 



. 3 



For the matter component we impose as in the previous case 8m izi) = tx — an( l nn d #a/(zi) = — ( 1+2a A — A ^ ( 1+z ) — j h 

L+Zi \ 1+2/3 A(i-A) (i+z)- -a J {L+Z *> 

The left panel of Fig. [5]shows the evolution of 8m, dot-dashed (red) lines, for different options for both Sx(zi) and 
9x(zi). Each different dot-dashed (red) line depicts two initial conditions, a — and /3 = and a = 1 and j3 = 10 
(they practically overlap). In this scenario, the final result does not depend on the chosen initial conditions. 
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3. Interacting version of model 2 



and 6i 



For the matter component we obtain as in the previous section lVI A l\ Si\r (zA = jq^yr cw±«-i u M — — \ n M0 +p i 
For the DE component we take different options for both 5x{zi) and 8x(zi). The dashed (green) lines of the right 
panel of Fig. [S] correspond to the noninteracting scenario. The solid (blue) line is for the ACDM model, shown for 
comparison purposes. Each different dashed (green) line depicts, from top to bottom in the right, a = and f3 — 0, 
a = 1 and (3 = 1 and a = 1 and /3 = 10. Notice that the interacting scenarios exhibit a dependence on the initial 
conditions chosen for the DE component. 



H 



4- <? version of model 2 
and 6 M (zi) = - 



For the matter component we find 8m{zi) = yqrj- an d 0Ai{zi) = — pjfrp • ^he right panel of Fig. \5\ shows 
the evolution of 5m, with dot-dashed (red) lines, for various options for both 8x(zi) and 8x(zi). Each different 
dot-dashed (red) line depicts the initial conditions, a — and ft — 0, and a — 1 and /? = 10, though they practically 
overlap. Again, the final result does not depend on the initial conditions chosen. 



B. Results 

Figure [5] shows the numerical solution for 5m for model 1 (left panel), and for model 2 (right panel), in both cases 
for k = 0.01 h Mpc -1 . However, the outcome is quite independent of the wavenumber k in the range 0.001 hMpc' 1 < 
k < 0.1 hMpc' 1 that includes all the interesting scales under consideration. In the interacting case, the matter 
density perturbations do not depend much on initial conditions imposed on the 5x, but they do on the initial 
conditions on 6x- Notice that in any case, the matter density perturbations clearly differ in both scenarios, the 
interacting and the c 2 one. The most favored (the closer to ACDM), at least in the two models studied here, are the 
c 2 scenarios, since low density perturbations at z ~ 10 can be problematic for the large scale structure formation. 
To confront it with observations, we resort to the growth function, / = e( In <5m /d hi a [3], and the observational 




Figure 5. Left panel: The evolution of DM density perturbations versus redshift in model 1. Right panel: The same for 
model 2. In plotting the graphs the dashed (green) lines describe the matter density perturbations of the interacting scenario, 
the dot-dashed (red) lines the c 2 , and the solid (blue) line ACDM for comparison purposes. Different initial conditions are 
used for the interacting versions; from top to bottom at the right, a = and /3 = 0, a = 1 and = 1 and a = 1 and /9 = 10. 
For the c 2 , two different sets of initial conditions are used, but the corresponding graphs practically overlap each other (a = 
and /9 = 0, and a = 1 and /3 = 10). Notice that in the right panel, the c 2 scenario overlaps ACDM, as it behaves as a toCDM 
model with w = 0.99. 



data borrowed from [47]. In Fig. HJ] we can see that for both models (left panel model 1, right panel model 2) the 
noninteracting version fits the data better, as said before, especially at the present epoch. This approach was just 
done to show that interacting versions and c 2 sharing the same H(z), and so indistinguishable at the background 
level, can evolve diversely at the perturbative level. However, to definitively discard, or validate, the interacting 
model, an exact first order perturbation together with the study of matter and radiation power spectrum appears 
necessary. This will be the subject of a future work. 
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Fi gUI6 6. Left panel: the evolution of the growth function, / = — ry — — , versus redshift for model 1. Right panel: the same 
for model 2. In plotting the graphs the dashed (green) lines describe the growth function of the interacting scenario, the 
dot-dashed (red) lines the c 2 , and the solid (blue) line ACDM for comparison purposes. Different initial conditions are used 
for the interacting versions; from top to bottom a = and = 0, a = 1 and /3 = 1 and a = 1 and /3 = 10. For the c 2 , two 
different sets of initial conditions are used, but they practically overlap each other (a = and j3 = 0, and a = 1 and (3 = 10). 
Notice that in the right panel, the c 2 scenario overlaps ACDM, as it behaves as a uiCDM model with w = 0.99. Observational 
data is borrowed from [471 ], 



VII. CONCLUSIONS 

As we have seen, at background level, holographic interacting models whose IR cutoff is set by the Hubble's length 
can be viewed as noninteracting ones whose holographic parameter c 2 is not constant but varies slowly. Because 
they share identical Hubble function, they are not distinguishable at the background level. However, this degeneracy 
can be broken at the perturbative level, since both energy components, i.e., DM and DE, evolve diversely in each 
scenario. The interacting model and the noninteracting one (section llVp fit the geometrical data rather well. At the 
perturbative level, the noninteracting scenario is favored by large scale structure formation. In the second model 
(Sec. [V}, which behaves as a spatially flat wiCDM model, both interpretations, interacting and c 2 (t), fit the data well 
and the values of the parameters seem reasonable. It contains the ACDM as a limiting case but with the coincidence 
problem alleviated in the interacting scenario. To discriminate between both interpretations, at the perturbative 
level, the noninteracting scenario appears favored. This should not be surprising, since it mimics rather well the 
ACDM model. 

In general, the noninteracting versions seem to be favored by the structure formation; however, the interacting cannot 
be ruled out just at this point. To go deeper in the matter, a full-fledged perturbative analysis should be undertaken 
and will be the subject of a future work. 

ACKNOWLEDGMENTS 

We are indebted to Diego Pavon for comments, advice and carefully reading the manuscript. We are also 
indebted to Gaetano Vilasi, Ninfa Radicella and Fernando Atrio-Barandela for fruitful discussions. I.D. and L.P. 
were supported by an INFN/MICINN collaboration under Grant. No. AIC10-D-000481. This work was partially 
supported by the "Ministerio Espafiol de Ciencia e innovacion" under Grant No. FIS2009-13370-C02-01, by the 
"Direccio de Recerca de la Generalitat de Catalunya" under Grant No. 2009SGR-00164 and by the Italian Ministero 
Istruzione Universit'a e Ricerca (MIUR) through the PRIN 2008 grant. I.D. was funded by the "Univesitat Autonoma 
de Barcelona" through a PIF fellowship. 



Appendix A: Hubble functions considering the radiation component 

To constrain the model with CMB data, it is necessary to take into account the radiation component, since 
we need to describe the Universe at the last scattering surface, z* ~ 1090, where the said component is no longer 
negligible. So, expressions like (|18|) and (I27p are not accurate at the last scattering epoch, much less so at earlier 
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times. The presence of radiation invalidates the expression pu — 3AfJ,(l — c 2 )H 2 , so to obtain the Hubble function 
we rewrite the second Friedmann equation as 

H 3 

■^2 = _ 5 ^ + w x^x +w R n R ) , (Al) 

where the subscript R stands for radiation, and the EoS parameter wx does not coincide with the w of Sec. II VI and 
IVl Differentiating Eq.® and using Eqs. (|3]2) and (|A1|I . we obtain 



Q wr ( H 



2 



w * = - (i-Z )PX H + ir{ir) ^+*) 4 - ( A2 ) 

With an expression for the EoS parameter of DE, Eq. (|A1|) has no analytical solution. However, we can consider 
two different integration regions, one for z < 50 and another for z > 50. In the former, the Universe is DM and DE 
dominated; in the latter only the higher z terms contribute. 

In the example considered in Sec. IIVI Q = Tpx, and in the first integration region, z < 50, the second term 
in Eq. (|A2f) contributes less than 1%, so the Hubble parameter is just as in Eq. (|18[) . an d the first Friedmann 
equation can be approximated by 

rr \ 2 

— ^A 2 + 2A(l-A)(l + z)i +(l-A) 2 (l + zf. (A3) 

H °J z<50 



JJS2 

H 



For the second region, z > 50, the first term in Eq. (|A2[) contributes no more than 1%, so the first Friedman equation 
is approximated by 

~ (l-A) 2 (l + z) 3 + 2fl R0 (l + z) i . (A4) 

z>50 

As the two first terms in Eq. (|A3[l are negligible in the region z > 50, and the second term in Eq. (|A4[l is trifling 
when z > 50, we can just consider the Hubble function, after redefining 2Q.ro — > Qro as 

H_ V I . . . ,i\ 2 

Ho 



A+(l-Ay(l + z)2) +Q R0 {l + z)\ (A5) 



This Hubble function is the same in both scenarios, the interacting and the c 2 , as when we ignored the radiation 
component. 

Proceeding as before, in the case of the model 2 (Sec. [Vj) we obtain 



Ho 



fi*ro(l + z) + «xo(l + 2 ) + n «o(l + z) ■ (A6) 



In both cases, Qro ~ f2 7 o(l + 0.0227^//) as described in [301 ] . where N e ff = 3.04 is the effective number of neutrino 
families. 



Appendix B: Lagrangian formulation of the models from sec. IIVI 
1. Interacting model with 1 scalar field 

In Sec. [T] we have seen that this kind of interacting models, can be described by just one fluid (as long as the 
evolution of the densities of DE and DM are identical), whose effective EoS parameter is 

with A a positive constant. The total effective pressure P — w e ffp is 

P = -3^ A M P H pi , (B2) 

where P and p = pm + px are the dark sector pressure and energy density from Eqs. (QJ and (|BI|) . For a standard 
scalar field cj> minimally coupled to gravity, the action is defined by 

S = j d^x^—g (j + L+fa X )) (B3) 
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where x = ^g^ud^^d" '<f> is the kinetic term. Assuming homogeneity and isotropy, the energy (Hamiltonian) density 
and the pressure (Lagrangian density) are 

P=\tf + V{<p) and P=l0 2 _y(0). (B4) 

From the last two equations the potential and the kinetic term are 

2 = p + P and V(<t>) = ^(p-P) , (B5) 

respectively. Using again Eqs.fT}, (|18[) and (|B2|) we obtain P(z) and with (|B5|l we obtain 

O / p2 \~5 

^- 1 +P) dP. (B6) 



^ 3,4^0 V3A 2 M 2 # 2 
After integration we get 

T n 3A 2 M 2 P H 2 ( A -iJ^frtt-M i.fjk^-MV 

L = P = ^ - Uf r'p + e V m p j , (B7) 

where we have defined (j>o as the value of the scalar field at the maximum pressure, i.e., when z — > — 1. It is interesting 
to notice that in the absence of interaction, A — 0, the pressure of the dark sector is that of standard matter. From 
(|B7ll and (|B5|) we find the potential 

3A 2 M 2 H 2 f A I^-WV L "l/? (Mo) nl -TtJ^»~M , 
1/(0) = ^ - I 4 + e V m p I 16e V a/ p + 24e V m p +1 

So we can describe the whole dark sector by a standard scalar field, that has an effective sound speed c 2 = ^ = 1, 

even if its adiabatic sound speed c 2 = is negative due to the interaction. Had we considered a general k-essence 
Lagrangian density [48|, the squared sound speed would have resulted negative. 



2. c 2 (t) model with 1 scalar field 

Here we obtain a Lagrangian density for the DE in the case of the holographic c 2 model. This is motivated to 
show that the two models of Sec. IIVI described at the background level by the same H(z), differ not only because 
the energy densities evolve differently but also because their Lagrangians are diverse. 
To begin with, the DE density can be expressed as a function of the pressure 

1 2 tt2 ti r2 



3A 2 H$Mp - 2P . (B8) 



Proceeding as before we obtain 



L = -3A'Mf.H$ sec 2 [ J ^ ^—f^ ) (B9) 



%A 2 M 2 P Hl ^ 2 / / 3 - 0o 



and 

^)-^^tan 2 (^^i] (BIO) 




As in the previous case, for a k-essence Lagrangian we would have obtained c s < 



[1] I. Duran, D. Pavon and W. Zimdahl, JCAP07, 018 (2010). 

[2] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116 . 1009 (1998); S. Perlmutter et al. [Supernova 

Cosmology Project Collaboration], Astrophys. J. 517 , 565 (1999). 
[3] T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis. [arXiv:1106.2476l [astro-ph.CO]. 

[4] J. A. Friemann, M.S. Turner and D. Huterer, Ann. Rev. Astron. Astrophys. 46, 385 (2008); R. Durrer and R. Maartens, 
Gen. Relativ. Grav. 40, 301; L. Amendola and Tsujikawa, Dark Energy. Theory and Observations (CUP, Cambridge, 
2010). 

[5] J. -H. He, B. Wang, E. Abdalla and D. Pavon, JCAP 1012 (2010) 022 



15 



[6] E. Abdalla, R. Abramao and J. C. C. Souza, Phys. Rev. D 82, 023508 (2010). 

[7] E. R. M. Tarrant, C. van de Bruck, E. J. Copeland and A. M. Gree n, Phys. Rev. D 85, 023503 (2012). 

[8] G. 't Hooft, "Dimensional reduction in quantum gravity", preprint gr-qc/9310026| L. Susskind, J. Math. Phys. (N.Y.) 
36, 6377 (1995). 

[9] A. G. Cohen, D.B. Kaplan and A.E. Nelson, Phys. Rev. Lett. 82, 4971 (1999); 

[10] S. D. H. Hsu, Phys. Lett. B 594, 13 (2004); D. Pavon and W. Zimdahl, Phys. Lett. B 628, 206 (2005). 

[11] M. Li, Phys. Lett. B 603, 1 (2004). 

[12] C. Gao, F. Wu, X. Chen, and Y.G. Shen, Phys. Rev. D 79, 043511 (2009). 

[13] L. Xu, W. Li, and J. Lu, Mod. Phys. Lett. A 24, 1355 (2009). 

[14] M. Suwa, T. Nihei, Phys. Rev. D 81, 023519 (2010). 

[15] I. Duran and D. Pavon, Phys. Rev. D 83, 023504 (2011). 

[16] S. Lepe and F. Pena, European Physical Journal (in the press), arXiv:10052.2180 [hepth]. 

[17] R. Brustein, "Cosmological entropy bounds", in String Theories and Fundamental Interactions, Lecture Notes Physics 

737 , 619 (2008), |hep-th/0702108| 

[18] W. Zimdahl and D. Pavon, Class. Quantum Grav. 24, 5461 (2007). 

[19] D. Pavon and B. Wang, Gen. Relativ. Grav. 41, 1 (2009). 

[20] S. del Campo, J. C. Fabris, R. Herrera and W. Zimdahl Phys. Rev. D 83, 123006 (2011). 

[21] D. Pavon, W. Zimdahl, Phys. Lett. B 628, 206 (2005) . 

[22] D. Pavon, J. Phys. A 40, 6865 (2007). 

[23] V. H. Cardenas, R. G. Perez, Class. Quant. Grav. 27, 235003 (2010). 

[24] B. Guberina, R. Horvat, H. Nikolic, JCAP07, 012 (2007). 

[25] L. Xu, JCAP09, 016 (2009). 

[26] H. Wei, Nucl. Phys. B 819, 210 (2009). 

[27] N. Radicella, D. Pavon, JCAP10, 005 (2010). 

[28] R. Amanullah et al. (The Supernova Cosmology Project), Astrophys. J. 716, 712 (2010). 

[29] W.J. Percival et al, Mon. Not. R. Astron. Soc. 401, 2148 (2010). 

[30] E. Komatsu et al, Astrophys. J. Supl. Series, 192, 18 (2011). 

[31] S.W. Allen et al, Mon. Not. R. Astron. Soc. 383, 879 (2008). 

[32] A.G. Riess, et al, Astrophys. J. 699, 539 (2009). 

[33] E. Gaztanaga, A. Cabre, and L. Hui, Mon. Not. R. Astron. Soc. 399, 1663 (2009). 

[34] J. Simon, L. Verde, and R. Jimenez, Phys. Rev. D 71, 123001 (2005). 

[35] D. Stern, R. Jimenez, L. Verde, M. Kamionkowski, and S.A. Stanford, JCAP02 (2010) 008. 

[36] S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 72, 123519 (2005). 

[37] D.J. Eisenstein et al. [DSS Collaboration], Astrophys. J. 633, 560 (2005). 

[38] E. A. Kazin et al, Astrophys. J. 710, 1444 (2010). 

[39] M. Kowalski et al, Astrophys. J. 686, 749 (2008). 

[40] J. Sollerman et al, Astrophys. J. 703, 1374 (2009). 

[41] S. M. D. White, J. F. Navarro, A. Evrard, and C. S. Frenk, Nature 366, 429 (1993). 

[42] G. Hasinger et al, Astrophys. J.573, L77 (2002); A.C.S. Friaga, J.S. Alcaniz, and J.A.S. Lima, Mon. Not. R. Astron. 

Soc. 362, 1295 (2005). 

[43] M.C. Bento, O. Bertolami, and A. A. Sen, Phys. Rev. D 66, 043507 (2002). 

[44] H. Kodama and M. Sasaki, Prog. Theor. Phys. Suppl. 78 1 (1984); 

[45] J. Valiviita, E. Majerotto and R. Maartens, JCAP 07 020 (2008); G. Caldera-Cabral, R. Maartens and B. M. Schaefer, 

JCAP 07 027 (2009). 

[46] L. Wang and P. J. Steinhardt, Astrophys. J.508, 483 (1998). 

[47] Y. Gong, Phys. Rev. D 78 123010 (2008). 

[48] C. Armendariz-Picon, V. Mukhanov and P. J. Steinhardt, Phys. Rev. D 63, 103510 (2001). 



